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Abstract 

Monopole-antimonopole pair (MAP) with both electric and magnetic 
charges are presented. The MAP possess opposite magnetic charges but 
they carry the same electric charges. These stationary MAP dyon solutions 
possess finite energy but they do not satisfy the first order Bogomol'nyi 
equations and are not BPS solutions. They are axially symmetric solutions 
and are characterized by a parameter, — 1 < r/ < 1 which determines the 
net electric charges of these MAP dyons. These dyon solutions are solved 
numerically when the magnetic charges of the dipoles are n = ±1, ±2 and 
when the strength of the Higgs field potential A = 0, 1. When A = 0, the 
time component of the gauge field potential is parallel to the Higgs field 
in isospin space and the MAP separation distance, total energy and net 
<^> electric charge increase exponentially fast to infinity when rj approaches ±1. 

T 7"! However when A = 1, all these three quantities approach a finite critical 

value as r\ approaches ±1. 

"X 

1 Introduction 

The SU(2) Yang-Mills-Higgs (YMH) field theory in 3 + 1 dimensions, with the 
Higgs field in the adjoint representation possess magnetic monopoles solutions 
[T]-[2] as well as dyons solutions |3j. The 't Hooft-Polyakov monopole solution 
with non zero Higgs mass and self-interaction is the first monopole solution that 
possess finite energy and this solution can acquire an electric charge to become 
a dyon [3J. These numerical, spherically symmetric monopole and dyon solutions 
of unit magnetic charge are invariant under a U(l) subgroup of the local SU(2) 
gauge group. Exact monopoles and multimonopoles solutions exist only in the 
Bogomornyi-Prasad-Sommerfield (BPS) limit [2]. Outside of this limit, when the 
Higgs field potential is non vanishing only numerical solutions exist. The exact 
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dyon solutions found by Prasad and Sommerfield [3] in the BPS limit is stable as 
they are the absolute minima of the energy [3]. 

The word, "dyon" was coined by J. Schwinger in 1969 j5] for a particle that 
possesses both magnetic and electrical charges. A dyon with a fixed magnetic 
charge can possesses varying electric charges. The dyons solutions of Julia and 
Zee j3] are time independent solutions that possess non vanishing kinetic energy. 
Dyons with axial symmetry were constructed by B. Hartmann et al. [6] when 
< r] < 1. These axial dyons are actually 't-Hooft-Polyakov monopoles that 
possess magnetic charges, n = 1, 2, 3 and non vanishing electric charges when the 
parameter rj is nonzero. It was found that when the strength of the Higgs potential 
A is non vanishing, the total energy and total electric charge of the system approach 
finite critical values when the parameter r] approaches 1. However when A = 0, 
the total energy and total electric charge of the system approach infinity when 
the parameter rj approaches 1. We found that the MAP dyon possesses similar 
characteristics. 

In this paper, we would like to present the axially symmetric MAP dyons by 
introducing electric charge to the system. The procedure of introducing electric 
charges to the monopole solution is standard as was first shown by Julia and Zee 
in 1975 [3]. The monopole and antimonopole carry the same electric charges and 
hence experience a repulsion force due to the electric charge. These MAP dyon 
solutions possess finite energy but they do not satisfy the first order Bogomol'nyi 
equations and are not BPS solutions. They are characterized by a parameter, 
— 1 < rj < 1 which determines the net electric charges of these MAP dyons. 
When 7] is positive (negative), the electric charges carried by both monopole and 
antimonopole are positive (negative). These dyons solutions are solved numerically 
when the magnetic charges of the dipoles are n — ±1, ±2 and when the strength 
of the Higgs field potential A = 0, 1. When A = 0, the time component of the 
gauge field potential is parallel to the Higgs field in isospin space and the MAP 
separation distance, total energy and net electric charge of the MAP dyons increase 
exponentially fast to infinity when 7] approaches ±1. However when A = 1, all these 
three quantities approach a critical finite value when t] approaches ±1. 

We briefly review the SU(2) Yang-Mills-Higgs field theory and the ansatz in the 
next section. We present our MAP dyon solutions and the results of our numerical 
calculations in section 3. We end with some comments in section 4. 

2 The SU(2) Yang-Mills-Higgs Theory 

The SU(2) YMH Lagrangian in 3+1 dimensions is given by 

£ = - l F"F a ^ - 1 D^ a D u ^ a - 1 A($ a $ a - — ) 2 , (1) 
4 ^ 2 4 A 

where fi is the Higgs field mass, A is the strength of the Higgs potential and 
£ = fi/y/X is the vacuum expectation value of the Higgs field. The Lagrangian 
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([!]) is gauge invariant under the set of independent local SU(2) transformations at 
each space-time point. The covariant derivative of the Higgs field and the gauge 
field strength tensor are given respectively by 

D^ a = d^ a + e abc A b ^ c , (2) 
F; v = d,A a v -d v Al + e abc AlAl. (3) 

Since the gauge field coupling constant, can be scaled away, we set it to one without 
any loss of generality. The metric used is g^v = ( — h ++)• The SU(2) internal 
group indices a, b, c run from 1 to 3 and the spatial indices are fi,u,a = 0, 1, 2, and 
3 in Minkowski space. 

The equations of motion that follow from the Lagrangian ([I]) are 

D»F* U = d^F* v + e abc A b »F° u = e abc ^ b D^ c , 

D»D^ a = A$ a ($ b $ fe - ^-). (4) 

A 

Non-BPS solutions to the YMH theory are obtained by solving the second order 
differential equations of motion Q, whereas BPS solutions can be more easily 
obtained by solving the Bogomol'nyi equations, 

Bf ± Di<S> a = 0, (5) 

which is of first order. 

The Abelian electromagnetic field tensor as proposed by 't Hooft [1] is 

= $ a F^-e abc $ a D^ b D u $ c 

= d^-d^-e^&d^dA , (6) 

where A^ = § a A a ^ <l a = $ a /|$|, |$| = y/$ a $ a . Hence the 't Hooft electric field 
is Ei = F i0 , and the 't Hooft magnetic field is Bi = —\eijkFj kl where the indices, 
i,j, k = 1,2,3. The topological magnetic current, which is also the topological 
current density of the system is [7j 

K = ^~ e, upa e abc d^ b d°&. (7) 

Therefore the corresponding conserved topological magnetic charge is 

M — J d 3 x k = — jd 2 (Ji B,. (8) 

In the BPS limit when the Higgs potential vanishes, the energy can be written 
in the form 

E = tJ d l (B^ a ) d 3 x + J -(B? ± D^ a f d 3 x 

= T J d^B^) d 3 x = 4vrM-^, (9) 

where M is the "topological charge" when the vacuum expectation value of the 
Higgs field, is non zero coupled with some non-trivial topological structure of 
the fields at large r. 
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3 The Dyons 



3.1 The Ansatz 

The time independent gauge fields and Higgs field that leads to the MAP dyon 
solutions are given respectively by, [8] 

At = - 1 Mr, $)ni>i + 6)n a e $i + l R 1 (r, e)h%fi - l R 2 (r, 9)n a r k 

ry r ry ry r ry 

A a = A 1 (r,9)n a r +A 2 (r,9)n a , 

$ a = § x {r,e)h a r + §2{r,e)h a e , (10) 

where the spatial spherical coordinate orthonormal unit vectors are 

fi = sin0 cos0 Sn + sin0 sin0 5 i2 + cos0 <5 i3 , 
§i = cos 9 cos <p 5n + cos 9 sin0 5 i2 — sin^ 5j 3 , 
4>i = - sin0 5n + cos0 5 i2 , (11) 

and the isospin coordinate orthonormal unit vectors are 

n a r = sin 9 cos ncj) 5^ ~\~ sin 9 sin ncp 5^ ~\~ cos 9 5^, 
•hg = cos 9 cos n<f) 5^ + cos 9 sin n(f) 5% — sin 9 5%, 
h a < = —sinn<j) 5^ + cosncf) where n > 1. (12) 



The constant parameter — 1 < r\ < 1 in Eq. (10) can also be written as r\ = tanh7. 
The 0-winding number n is a natural number. The MAP solutions exist only when 
n = 1,2. When n = 3, the MAP disappeared and a vortex ring is formed instead 
[Sj. Hence in our numerical calculation, we solved for the MAP dyon solutions 
only when n = 1, 2. 

3.2 The Solutions 



In order to solve the for the MAP dyons solutions, the ansatz (10) is substituted 
into the equations of motion Q and the 15 equations of motion are reduced to 
eight coupled second order partial differential equations. These eight equations are 
solved asymptotically first at small and then at large distances. In these asymptotic 
regions, the time component of the gauge field and the Higgs field are parallel in 
the isospin space, that is $i = A\ and $ 2 = A 2 - The asymptotic solutions and 
boundary conditions at small distances that will give rise to finite energy MAP 
dyons solutions are 

Vu(O,0) = Ra{0,0) = 0, A = 1,2; 
sin0 $i(O,0) + cos0 $ 2 (O,0) = 
sin0 A(O,0) + cos0 A 2 (O,0) = 
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0_ 
dr 

— i 
dr 



(cos9 $i(r, 9) - sm.9 $ 2 (r, 9)) 
cos 9 Ai{r,9)-Bin9 A 2 (r,9)) 



r=0 



r=0 




0. 



At large distances, the asymptotic solution can be written as 



(13) 



ipx(oo, 9) = 2, ^2(00, 9) = 2n, Ra(oc, 
$1(00, 9) = £cos#, $2(00, 
Ai(oo, 9) = r)£ cos 9, ^.2(00, 



0, 

£sin6>, 

^sin6» (14) 



Ra(t, 
d e tlJA(r,9)\e- 



?)|«->o, 7r = $ 2 (r,i 
,0, 7r = <%$i(r,#) 



|e->o, TT = ^4-2(^, 6 l )|e^o, tt 
->o, 7r = <9e^i(r, 6 l )| e ^ , 7r 



0. 
0. 



(15) 



where A = 1,2, and the expectation value £ = ^% — 1. The numerical MAP dyons 
solutions connecting the asymptotic solutions (13) to (14) and subjected to the 

(16) 



boundary conditions (15) together with the gauge fixing condition [8] 

d r Ri — rdgipi = 0, 



were solved using the Maple 12 and MatLab R2009a softwares. The second order 
equations of motion Q which are reduced to eight partial differential equations 
with the ansatz (10) are then transformed into a system of nonlinear equations 



using the finite difference approximation. This system of nonlinear equations are 
then discretized on a non-equidistant grid of size 70 x 60 covering the integration 
regions < x < 1 and < 9 < 71. Here x is the finite interval compactified 
coordinate given by x = The partial derivative with respect to the radial 



coordinate is then replaced accordingly by d r — > (1- 



2(1 -x 



,3_9 

dx 



-x) 2 d s and — > 



a 2 



First of all, we used Maple to find the Jacobian sparsity pattern for 



the system of nonlinear equations. After that we provide this information to 
Matlab to run the numerical computation. The system of nonlinear equations are 
then solved numerically using the trust-region-reflective algorithm by providing 
the solver with good initial guess. 

The second order equations of motion Eq. Q were solved when the 0-winding 
number n = 1, 2 and with zero Higgs potential, that is A = fi = but with nonzero 
expectation value £ = 1 and finally with A = // = 1 and expectation value £ = 1. 

In the Prasad-Sommerfield limit where A = 0, the ansatz (10) can be simplified 
by letting 

Ax{r, 9) = 77$ x (r, 9) and A 2 (r, 9) = ^ 2 (r, 9) (17) 
hence making the gauge potential Aq parallel to the Higgs field $ a in isospin space. 
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Upon substituting ansatz (10) and (17) into the equations of motion Q, the 
resulting equations of motion are simplified to, 



D l F a 

13 



Qi F a + e abc A U F c 



D l D, } <$> a = 0, 



(18) 



where $ a = yjl — rf $ a = sech7 $ a , the covariant time derivative D^ a = 0, and 
the electric field E? = F® = —r]Di^ a . The time component of the 't Hooft gauge 
potential becomes, A = Hence the 't Hooft electric field is non zero when 

7] 7^ and the monopole becomes a dyon. 



3.3 The Magnetic Field, Magnetic Charge, and MAP Sep- 
aration 
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4.1970 


1.6948 





1.7944 


2.9514 


0.20 


4.3573 


4.2840 


1.7298 


7.6013 


1.8312 


3.0126 


0.40 


9.3156 


4.5806 


1.8496 


16.2519 


1.9588 


3.2216 


0.60 


16.0064 


5.2502 


2.1196 


27.9272 


2.2476 


3.6928 


0.80 


28.4471 


7.0118 


2.8268 


49.6461 


3.0084 


4.9274 


0.90 


44.0378 


9.6884 


3.8908 


76.9414 


4.1712 


6.7909 


0.97 


84.8568 


17.6718 


6.9571 


150.4439 


7.7512 


12.3184 


0.98 


105.2895 


21.9170 


8.5435 


183.2625 


9.6248 


14.8613 



Table 1: Values of Q(n, X,rj), d(n,X,r]), and E(n,X,r]), when A = and n = 1,2. 
When 77 — y 1, Q, d, and E diverges. 

The 't Hooft magnetic field Bi is made up of two parts, the gauge part and the 
Higgs part, 

B, = Bf + B? , 

Bf = -e ijk d 3 A k and B* = ^e ijh (? bc $ a d j $ b d k $ c . (19) 

With some calculations, the gauge part of the magnetic field is given by 

Bf = -ne ijk dj sin k d k (f>, 

sm9 ( $ 2 $1 \ . , 

where, shift = — — \ ip 2 j^y - R2 t^t I , (20) 

To calculate for the 't Hooft magnetic field Bf , we rewrite the Higgs field from 
the spherical to the Cartesian coordinate system, 

$ a = $1 n a r + $ 2 h% + $ 3 h% 

= $! 5 al + $ 2 5 a2 + $ 3 <5 a3 , (21) 
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A = 1 


n = 1 


n = 2 
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E 
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3.2572 


2.3635 





1.4958 


4.8582 


0.20 


2.7387 


3.2956 


2.3854 


3.7496 


1.5122 


4.8882 


0.40 


5.6103 


3.4204 


2.4543 


7.5969 


1.5598 


4.9805 


0.60 


8.7866 


3.6652 


2.5813 


11.6555 


1.6518 


5.1428 


0.80 


12.5557 


4.1264 


2.7924 


16.0816 


1.8068 


5.3906 


0.90 


14.8390 


4.5218 


2.9470 


18.5080 


1.9210 


5.5551 


0.97 


16.7087 


4.9368 


3.0860 


20.3277 


2.0224 


5.6908 


0.98 


17.0015 


5.0116 


3.1087 


20.5977 


2.0384 


5.7118 


1.00 


17.6100 


5.1760 


3.1570 


21.1500 


2.0720 


5.7550 



Table 2: Values of Q(n, X, rj), d(n,\,r)), and E(n,X,T]), when A = 1 and n = 1,2. 
When T] — > 1, Q, d, and E tends to their respective maximum critical values. 

where $i = sin#cos7i(/> $i + cos#cosn0 $2 — sinn0 $3 = |$| cos a sin /3 
$2 — sin#sinn0 $ x + cos#sinn0 $2 + cosn0 <3> 3 = |$| cosacos/? 
$3 = cos# $1 - sin# $ 2 = |$| sina. (22) 

The Higgs unit vector can be simplified to 

<f> a = cos a sin /3 5 al + cos a cos /3 5 a2 + sin a 5 a3 , 

sin a = t— cos 6* — -— sm9, (3 = neb, (23) 

$ $ 2 

and the 't Hooft magnetic field is reduced to only the fj and Q{ components, 

B? = -ne ijk d j sin ad k (f) 

n [ d sin a 1 n [ d sin a 1 - . . 

r l sm t) [ o(J ) r sin # [ or J 

From our numerical calulation, we managed to plot magnetic field lines of the 
MAP dyons when n = 1 for values of A = 0, 1 and i] = 0.20, 0.90, 0.98 as in Figure 
[T] and when n = 2 for values of A = 0, 1 and i] = 0.20, 0.90, 0.98 as in Figure |2j 

From Eq.(|8]), the net magnetic charges enclosed by the upper and lower hemi- 
sphere at infinity are calculated to be, 



n . 

M, = sin a 

+ 2 

n . 

M_ = sm a 

2 



= +n, 

0,r— >oo 



-n, n = l,2, (25) 



respectively. Hence the magnetic charge in the upper and lower hemi-sphere is 
+n and — n respectively and the net magnetic charge of the MAP dyon is zero. 
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For the MAP dyon to exist, n = 1, 2 and when n = 3, the MAP dyon cease to 
exist and a vortex ring dyon is formed instead [5]. Following 't Hooft's definition 
for the electromagnetic field, Eq. ([6]), the magnetic charges of n = ±1,±2, are 
all located at the zeros of the Higgs field and the MAP separations d(n, A, 77) for 
different values of n, A, and 77 are as given in Table [T] and [2} 

From Figure [T] and Figure [2j we can generally conclude that when A = 1, the 
MAP separations are smaller than when A = for given values of rj and n as the 
interaction between Higgs field is attractive. When n = 2, the MAP separations 
are smaller than when n = 1 for fixed values of A and 77 as the MAP possess unlike 
magnetic charges that attract. When the value of 77 increases the electric charges 
also increases and for fixed values of A and n, the MAP separations get wider as 
the MAP possess similiar electric charges that repel. Hence the MAP separation 
d(n, A, 77) varies with n, A, and 77. From our numerical calculations, we found that 
when A = 0, d(n, 0, 77) increases with 77 from d(l, 0, 0) = 4.20 and d(2, 0, 0) = 1.79 
to infinity exponentially fast as 77 approaches one. However when A = 1, d(n, 1,77) 
increases with rj from d(l, 1, 0) = 3.26 and d(2, 1, 0) = 1.50 to finite critical values 
of d(l, 1, 1) = 5.18 and d(2, 1, 1) = 2.07 respectively as 77 approaches one. See 
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Magnetic Field t\ = 0.2. ), = 0, n = 2 



Electric Field. i] = 0.2, 1 = 0, n = 2 



Magnetic Field, n = 0.2,*.= 1,n = 2 
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Figure 2: The magnetic field lines and the electric field plots of the MAP dyons 
when n = 2 for (a) A = and (b) A = 1 when 77 = 0.20, 0.90, 0.98. 

Table Q] and 

3.4 The Electric Field and Electric Charge 

At spatial infinity in the Higgs vacuum, the electromagnetic field can be repre- 
sented uniquely by 



flU 



T 



(26) 



which corresponds to the unbroken U(l) symmetry. The corresponding electric 
and magnetic fields are then = J 7 ^ and Bi = — ^e^Tju- However according to 
Coleman [9], there is no unique way of representing the electromagnetic field in 
the region of the monopole outside the Higgs vacuum. One proposal was as given 
by 't Hooft as in Eq. (|6| which will reduced to Eq. (26) at spatial infinity in the 
Higgs vacuum. The 't Hooft's magnetic field has the special property, d l Bi = 
when |$| 7^ 0. Hence with 't Hooft's definition of the electromagnetic field, the 
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magnetic charges can only reside at the zeros of the Higgs field and the magnetic 
charges of the MAP dyon configurations are concentrated at two points along the 
z-axis where the Higgs field vanishes. Thus there are points singularities in 't 
Hooft's definition of Eq. (|6]), which makes it not possible to calculate and plot the 
magnetic and electric charge density distributions numerically with this definition. 

The 't Hooft electric field of the MAP dyon, E { = diA = ft + A 2 ^} 

corresponds to Eq. ( |26| ) in the Higgs vacuum. 2D vector plots of the electric field Ei 
are given in Figure [^when n = 1 and Figure [2] when n = 2, for r] = 0.20, 0.90, 0.98 
and (a) A = 0, (b) A = 1. Unlike the magnetic field, the electric field varies 
proportionally with the constant — 1 < r\ < 1. Hence the electric field can be 
switched off by setting rj = 0. As the parameter r] increases from —1 to zero, 
the electric field points radially inwards and when < r\ < 1, the electric field 
is radially outwards. Hence the electric charges of the both the monopole and 
antimonopole are positive (negative) when r] is positive (negative). 

Another proposal for the electromagnetic field, JF^ = F£ u {~}; wa s given by 



Bogolmolny [10J and Faddeev [UJ. Using this definition, Eq. (26), the magnetic 
charge density g = d l Bi and electric charge density q = d l £i distributions can be 
calculated and plotted numerically as shown in Figure [3] (a), (b) to Figure [6] (a), 
(b) when A = 0, 1 and n — 1, 2 respectively. The electric charge density, q = d l Si, 
of the MAP dyon solutions is solely positive (negative) throughout space when r] 
is positive (negative). When n = 1, the electric charge distribution of the MAP 
dyon is two spheres centered at the two point zeros of the Higgs field. Hence the 
soliton has a dumbbell shape structure. 

However when n = 2, the electric charge distribution is two horizontal toruses 
around the two zeros of the Higgs field as center with the z-axis as the symmetry 
axis. Contour plots of the electric charge density distribution when n = 1,2, 
A = 0, 1, and r] = 0.98 are given in Figure [3] (b) to Figure [6] (b). Contour plots 
of the magnetic charge density distribution, the Higgs field modulus, and energy 
density when n = 1, 2, A = 0, 1, and r] = 0.98 are also given in Figure [3] (a), (c), 
(d) to Figure M (a), (c), (d) respectively. 



From Gauss' law, the total electric charge of the dyon is given by, 

Q(n,\,r})= f Sih r 2 sin9 ded<p. (27) 

The different values of Q(n, X,rj) calculated for various values of n, A, and rj are 
as tabulated in Table [l] and [2j For fixed value of n — 1, 2 and A = 0, Q(r/) 
increases with rj from Q(n, 0, 0) = and approaches infinity exponentially fast as 
7] approaches one. However for fixed value of n = 1,2 and A = 1, Q(r]) increases 
with rj from Q(n, 1, 0) = to a maximum critical value of Q(l, 1, rj) = 17.61 and 
Q(2,l, rj) = 21.15 as rj approaches one. The four parametric graphs for MAP 
separation d(n, A, rj) versus total electric charge Q(n, A, rj) for n — 1,2 and A = 0, 1 
were plotted for —0.98 < r] < 0.98 when A = 0, Figure [7] (a) and for — 1 < r] < 1 
when A = 1, Figure [7] (b). The end points for the graphs when r] — > ±1 in Figure 
[7] (b) are (5.18, ±17.61) for n = 1 and (2.07, ±21.15) for n = 2 whereas there is 
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no end point for the graphs in Figure [7j (a) as d(n, A, 77) and Q(n, A, rj) diverge as 
77 — > ±1. The relationship between the MAP separation d(n,X,r]) and the total 
charge Q(n, A, 77) when A = 1 is given by 

d « 4.4360 x 10" 6 Q 4 - 9.2560 x 10~ 20 Q 3 + 0.4771 x 10~ 2 Q 2 

+ 2.7050 x 10~ 18 Q + 3.2630, n = 1, 

« 4.5710 x 10~ 7 Q 4 + 6.6840 x 10~ 20 Q 3 + 0.1084 x 10~ 2 Q 2 

- 2.7650 x 10~ 17 Q + 1.4960, n = 2. (28) 

3.5 The Energy 

Since the MAP dyon is not a BPS solution, its energy cannot be given by Eq. (J9]) 
but is given by 

1 f A 
E= - J {BfB? + E?E? + D^ a Di^ a + D $ a D $ a + -($ a $ a - £ 2 ) 2 } d 3 x. (29) 

The different values of E(n, X,1]) calculated when £ = 1 for different values of n, 
A, and 77 are as tabulated in Table [T] and [2] Similar to the MAP separations d 
and total electric charge Q, the total energy E approaches a finite critical value 
when 77 — > ±1 for A = 1 as shown in Figure [8] (b). When A = 0, the total energy 
E diverges exponentially fast as 77 — > ±1. The critical values of the total energy 
are £(1,1, ±1) = 3.157 and E(2, 1, ±1) = 5.755, see Figure @ (b). The energy 
density distribution is plotted 3D for n = 2, A = 1 in Figure [7] (a) when 77 = and 
(b) when 77 = 0.98. These plots show that the energy density of the MAP has a 
toroidal distribution around the two zeros of the Higgs field. 



4 Comments 



The ansatz (10) used to construct the dyon solutions is introduced in the same 
standard way as in the construction of the Julia- Zee dyons years ago [3]. Hence 
when the parameter 77 is zero, the electric field is switched off and only the magnetic 
field remains and the dyon becomes a monopole. 

The MAP separation d, total electric charge Q, and the total energy E of 
these MAP dyon solutions are infinite as 77 — > ±1 when A = or when the Higgs 
field $ a is parallel to the time component of the gauge field Aq in isospin space. 
However when $ a is not parallel to A$, that is when A = 1, all the three quantities, 
d(n,X,T]), Q(n,X,rj), and E(n,X,T]) approach critical values when 77 — > ±1 . The 
relationship between d, Q, and E are been studied and {d(n, 1, 77) — d(n, 1, 0)} is 
approximately proportional to Q 2 (n, 1, 77). 

These MAP dyon solutions are axially symmetric about the z-axis. Hence 
the MAP dyon is a magnetic dipole with both poles carrying the same electric 
charges but opposite magnetic charges. When n = 1, the electric and magnetic 
charges are two spheres with zeros of the Higgs field as center. The sphere in the 
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upper hemishphere possesses positive magnetic and electric charges whereas the 
sphere in the lower hemisphere possesses negative magnetic charge and positive 
electric charge when 77 is positive, Figure [3] and Figure |4| The MAP dyons possess 
dumbbell shape mass distribution. When r\ is negative, only the electric charges 
will reverse its sign. However when n — 2, the magnetic and electric charges are 
two toruses with the zeros of the Higgs field as their center. In this case, the two 
zeros of the Higgs field are multiple zeros as shown in Figure [5] (c) and Figure [6] (c). 
Hence the soliton has the unique energy density distribution as given in Figure [5] 
(d) and Figure [6] (d) when A = and A = 1 respectively. 

Our conclusion is that the MAP dyon when n = 2 is not a two points monopole- 
antimonopole but rather, it is composed of two toruses of the same electric charges 
but opposite magnetic charges, whereas when n — 1, the solition is composed of 
two spheres. 
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Magnetic Charge Density, r\ = 0.98, A = 0. n = 1 Electric Charge Density, n = 0.98, A = 0, n = 1 




(C) (d) 



Figure 3: Contour plots of (a) the magnetic charge density, (b) the electric charge 
density, (c) the Higgs field modulus, and (d) the energy density, when n — 1 for 
A = and rj = 0.98. 
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Magnetic Charge Density, n = 0.98, X = 1, n = 1 Electric Charge Density, n = 0-98, X = 1, n = 1 




Figure 4: Contour plots of (a) the magnetic charge density, (b) the electric charge 
density, (c) the Higgs field modulus, and (d) the energy density, when n — 1 for 
A = 1 and rj = 0.98. 
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Magnetic Charge Density, i] = 0.98, k = 0, n = 2 



Electric Charge Density, n = 0.98, k = 0. n = 




Figure 5: Contour plots of (a) the magnetic charge density, (b) the electric charge 
density, (c) the Higgs field modulus, and (d) the energy density, when n = 2 for 
A = and rj = 0.98. 
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Magnetic Charge Density, n = 0.98, K = 1, n = 2 Electric Charge Density, n = 0.98. = 1, n = 2 




(c) (d) 



Figure 6: Contour plots of (a) the magnetic charge density, (b) the electric charge 
density, (c) the Higgs field modulus, and (d) the energy density, when n = 2 for 
A = 1 and rj = 0.98. 
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Plot of electric charge versus distance between two dyons, n =1 , X = Plot of distance between two dyons versus electric charge, n = 1 , 1 = 1 
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Figure 7: Parametric plots of (a) d(n, 0, if) versus Q(n, 0, if) and (b) d(n, l,rj) versus 
Q(n, 1, if) when n — 1, 2 for — 1 < 77 < 1. 
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Energy Density, r\ = 0,\ = '\,n = 2 Energy Density, r\ = 0.98, X = 1 , n = 2 




(a) (b) 



Figure 9: 3D plots of the energy density distribution for n = 2, A = 1 when (a) 
r] = and (b) r] = 0.98. 
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